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O ■ ABSTRACT 

O ■ 

(N : 

pH . Accurate measurement of pulsar distances via astrometry using very long 

^ . baseline interferometry enables the improvement of Galactic electron density dis- 

^ . tribution models, improving distance estimates for the vast majority of pulsars 

CN . for which parallax measurements are unavailable. However, pulsars at southern 

I— i. declinations have been under-represented in previous interferometric astrometry 
campaigns, due to the dominance of northern hemisphere instrumentation for 

^ , astrometry. In order to redress this imbalance, we have conducted a two-year 

P_i! astrometric campaign targeting eight southern pulsars with the Australian Long 

Q ! Baseline Array. The program summarized in this paper has resulted in the mea- 

^ ! surement of seven new pulsar parallaxes, with success on objects down to a mean 

^ ! flux density of 800 /iJy at 1600 MHz. Our results highlight the substantial un- 

! certainties that remain when utilizing free electron density models for individual 

^ ! pulsar distances. Until this study, PSR J0630-2834 was beheved to convert 16% 

! of its spin-down energy into x-rays, but our measured parallax distance of 3321|q 

^ ! pc has revised this value to <1%. In contrast, PSR J0108-1431 was found to 

^ \ be almost a factor of two more distant than previously thought, making its con- 

^ \ version of spin-down energy to x-rays the most efficient known (>1%). The 8.5 

G^ ! second radio pulsar J2 144-3933 was found to be closer than previously predicted, 

^ \ making its apparent 1400 MHz radio luminosity the lowest of any known pulsar 

\ (20 /iJy kpc^). We have examined the growing population of neutron stars with 

^ \ accurate parallaxes to determine the effect of distance errors on the underlying 

5t ! neutron star velocity distribution, and find that typical distance errors may be 
biasing the estimated mean pulsar velocity upwards by 5%, and are likely to 
exaggerate the distribution's high-velocity tail. 

Subject headings: Techniques: interferometric — astrometry — pulsars: general 



^Centre for Astrophysics and Supercomputing, Swinburne University of Technology, Mail H39, P.O. Box 
218, Hawthorn, VIC 3122, Austraha 

^Department of Imaging and Applied Physics, Curtin University of Technology, Bentley, WA, Australia 

■^Australia Telescope National Facility, Epping, NSW, Australia 



- 2 - 



Introduction 



Studies of pulsars are often plagued by tlie uncertainty in distance-dependent quantities 
introduced by the reliance on dispersion measure (DM) distance estimates. Due to the 



non-uniform distribution of ionized material in the ISM on sub-kpc scales ( ICordes et al. 
20071 ). the correspondence between a pulsar's DM and distance is often uncertain. The 



two most widely used G alactic electron density distributi on models are those pr esented by 
Taylor fc Corded (119931 ). hereafter referred to as TC93, and lCordes fc Lazid (120021 ). hereafter 
referred to as NE2001. Although distances estimated from these models are typically quoted 
as being accurate to 20%, previous astrometric pulsar observations have shown that much 
greater errors ar e possible for individu al objects (e.g. a factor of 2.5 error in the distance for 
PSR B0656+14: ferisken et al.ll2003bh . 



Independent distance measures to pulsars are extremely valuable, both to enable the 
confident estimation of distance-dependent parameters for individual pulsars and to refine 
DM-based distance models for the remainde r of the pulsar popula tion. They can be pro- 
vided from HI absorption measurements (e.g . iKoribalski et al.lll995l ). association with other 



astrophysical objects (e.g. ICamilo et al 



20061 ). and ann ual geometric paralla x measurements 



made via timing (e.g. iHotan et al.l 120061 ) or VLBI (e.g. iBrisken et al.l 120021 ). Both scientific 



justifications are more acute at southern declinations, wher e prior to this observing program 



only two published VLBI pulsar parallaxes were available (IBailes et al.lll990l : iDodson et al. 



20031 ). despite the fact that 70% of known pulsars are located at a declination < 0°. Re- 
cently, a fu rther four parallaxes ha ve been published for pulsars with declinations between ° 
and —20° (IChatterjee et al.ll2009l ). At the present time, however. Galactic electron models 
still have fewer constraints in the south than at northern declinations, and even the most 
scientifically desirable targets such as the double pulsar J0737-3039A/B generally have only 
a DM distance for deriving luminosities, velocities, population sizes and so on. 

The eight pulsars targeted by this survey are described in Table [H Two of these (PSR 
J1559-4438 and PSR J2048-1616) were chosen simply as relatively bright nearby pulsars on 
which our data reduction techniques could be tested in the high signal-to-noise regime. The 
remainder were each chosen because an independent distance measure offered a resolution to 
an outstanding question surrounding the pulsar. Further constraints on target selection were 
imposed by the need for relatively uniform distribution around a range of right ascensions, 
due to observational logistics. Thus, the target selection in no way comprises an unbiased 
sample of southern hemisphere pulsars; however, as noted above, the paucity of existing 
parallax measurements in the south means that any additions will be extremely useful in 
constraining electron density distribution models at these declinations. 



The six pulsars selected on scientific merit can be broadly grouped into two categories. 
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The first encompasses those possessing unusual luminosity characteristics (low radio lumi- 
nosity pulsars PSR J0108-1431 and PSR J2144-3933, and the high x-ray luminosity pulsar 
PSR J0630-2834) where the primary motivation is to obtain a confident estimate of the 
pulsar luminosity. The second group encompasses binary pulsars used for tests of general 
relativity (GR; PSR J0437-4715, PSR J0737-3039A/B and PSR J2145-0750), where ac- 
curate distances and velocities are required to calc ulate kinematic contributions to pulsar 
timing terms (the Shklovskii effect; IShklovskiil 119701 ) . 



In this paper, we summarize our observational strategy in ^ and present the results of 
the survey in ^ Interpretation of the revised distance measure for each pulsar is presented 
in SI While the results for PSR J0437-4715 PSR J0737-3039A/B, and PSR J1559-4438 
have been presented previously ( Deller et al. 2008 . 2009al lbl). a summary is included here 
for completeness. The accuracy of TC93 and NE2001 predictions for pulsar distances are 
examined in SI along with an investigation into the likely impacts of distance model errors 
on derived pulsar population parameters. Our conclusions are presented in |j5l 



Observations 



The observational methodology employed has been discussed in detail in iDeller et al. 



( 2009b ). and a brief overview is presented here. Seven observing session^ spanned an 1 



month period between August 2006 and February 2008, with a subset of the eight pulsars 
observed at each epoch, depending on which were closest to parallax extrema. All observa- 
tions were made at 1650 MHz, with the exception of PSR J0437-4715 which was observed 
separately at 8400 MHz due to its brightness and narrow pulse profile. A phase reference 
cycle time of six minutes, allocated equally between calibrator and target, was used for all 
observations. 

Where possible, all six stations of the Australian Long Baseline Array (LBA) were used 
for observing. However, the Ceduna telescope does not possess a 1600 MHz receiver, and 
the use of a NASA DSN antennjfl in Tidbinbilla was possible only when telesc ope time was 



available, which was for roughly half of the observational program. Figure 1 of iDeller et al. 



( l2009bl ) shows representative uv coverage for 1650 MHz and 8400 MHz observations. While 
higher data rates were recorded at some stations, the final astrometric processing of the 
observations utilized four single polarization 16 MHz bands sampled at two bit precision. 



^As noted in IDeller et al.l (j2009bl ). it was necessary to discard the first of eight observed epochs (MJD 



53870) due to the probability of uncorrectable ionospheric model errors 
^The TOm and 34m antennas were used in different epochs 
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yielding a total data rate of 256 Mbps per telescope. 

The data were correlated using matched filteri ng (a more sophis ticated pulsar 'gate') 
on pulse profiles with the DiFX software correlator (jPeller et al.l 120071 ). which increases the 
signal to noise ratio of the correlated dataset by weighting the visibility data according to 
the predicted pulsar flux. The result ant correlated datase ts were processed with a pipeline 
utilising the Parseltong ue interface ( Kettenis et al.l bood ) to AIP^. This pipeline, which 
is described in detail in iDeller et al.l (l2009bl ). covered calibration, flagging, geometric and 
ionospheric model correction, and fringe-fitting, and corrected for reference source structure 
and pulsar scintillation. Figures [1] and [2] show the composite reference images used for the 
each phase reference source, and the position centroid used for each calibrator is given in 
Table [21 Station positions were taken from the global VLBI solution 2008b_astrcG, which 
in cludes new po si tion d eterminations for the Ceduna, Mopra and ATCA stations as described 



m 



Petrov et al.l (120091 ). The calibrated pulsar images produced at each epoch were then 



fit for position, and the multi-epoch position dataset thus generated was fit for reference 
position, proper motion and parallax. The astrometric fit at the final stage of the pipeline 
estimated and attempted to account for the presence of systematic errors in the pulsar 
position measurements. 



3. Results 

Of the eight pulsars observed, parallaxes were obtained for seven. The remaining pulsar 
PSR J2145-0750, appeared to suffer heavily from refractive scintillation, as was reported 



durin g a previous unsuccessful attempt at VLBI astrometry using the VLBA (IBrisken et al. 



20021 ) . In our observations, it was detected in only two usable epochs, and so while a proper 
motion can be derived by holding the parallax value fixed over a reasonable range of values, 
the results should be treated with some caution. 

Table [3] summarizes the results for the eight targe t pulsars. A l l error s are la values 



obtained using the "inclusive" fit approach described in iDeller et al.l (l2009bl ). which models 
the presence of systematic errors. The resultant reduced chi-squared value is 1.0 for all 
pulsars except PSR J0737-3039A/B (for which the initial fit gave a reduced chi-squared value 
of 0.79, and hence no additional systematic error component was required). All positions are 
specified at MJD 54100 (31 December 2006) and errors on positions are fit errors which do 



■^http:/ /www. aoc.nrao.edu/aips 

''http://lacerta.gsfc. nasa.gov/vlbi/solutions/2008b^stro/ 
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not include the several mas uncertainty typical to calibrator positions and constant phase 
referencing offsets. Table [3] lists the covariances between proper motion and parallax for 
each pulsar, and shows that the addition of future epochs could considerably improve the 
parallax precision for PSR J0437-4715 and PSR J2048-1616 in particular. 

Figures [314T0] plot the results for each pulsar, showing the astrometric fit overlaid on 
the measured positions with la error bars. In each figure, the pulsar position is shown in 
right ascension and declination in the top panel, while the middle and lower panel show 
right ascension plotted against time and declination plotted against time respectively. In the 
lower two panels, the best-fit proper motion has been subtracted for clarity. 



As discussed in lDeller et al.l ( l2009bl ). the optimal weighting for visibility data depends on 



the strength of the pulsar and th e magnitude of systematic errors. Following the guidelines 



developed in iDeller et al.l (l2009bl ). we used equally-weighted visibility data when the target 
could be detected in a single epoch using this weighting with a signal to noise ratio (SNR) 
greater than ~10; we have used sensitivity-weighted visibilities when this criterion could 
not be met. For the strongest pulsars, we have used uniform weighting when imaging to 
maximise resolution, but natural weighting was used for weaker pulsars where the use of 
uniform weighting has an adverse effect on single-epoch SNR. These choices are listed in 
Table O 



Along with the recent results of lChatterjee et al.l (120091 ) . our detection of seven southern 



hemisphere pulsar parallaxes increases the total number of VLBI pulsar parallaxes in the 
south to 12, and the total number of published VLBI pulsar parallaxes to 35. Based on 
the results presented here, it is possible to estimate the number of southern pulsars which 
remain accessible to the LBA for future parallax studies. The parallax error of 140 /xas 
obtained for PSR J0737-3039A/B, which possesses an ungated 1600 MHz flux of only 2-3 
mJy, shows that LBA parallax measurements for moderately faint pulsars are possible to 
distances of several kpc. A search of the ATNF pulsar catalogu^ reveals that there are 
over 30 suitable pulsars (declination < 0, 1400 MHz flux > 2.5 mJy and DM distance < 2 
kpc) that have not previously been the subject of VLBI astrometry. Future improvements 
to VLBI instrumentation should rapidly make even more pulsars accessible. Thus, future 
parallax observations with the LBA could and should continue to redress the imbalance of 
measured pulsar parallaxes at southern declinations. 

The implications of the measured parallax and proper motion for each pulsar are dis- 
cussed below in ^ along with the implications of the inferred error distribution for distance 
estimates based upon Galactic electron models and pulsar DMs. 



■''http://www.atnf.csiro.au/research/pulsar/psrcat/ 
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Discussion 



4.1. PSR J0108-1431 



PSR J0108-1431 was first reported by lXauris et al. ( 1994 ). and was postulated to be the 
nearest observed radio pulsar - its DM of 2.38 remains the lowest of any known pulsar. In 
the TC93 Galactic electron mo del, PSR J0108-1431 w as predicted to lie 130 pc from Earth 
- in the newer NE2001 model (jCordes fc Lazidl2002l ). the distance was revised to 184 pc, 
which if correct would double the pulsar luminosity, to a value of 300 yuJy kpc^ at 400 MHz 
(the lowest known value is that of PSR J0006+1834, which at a DM-derived distance of 700 
pc has a predicted 400 MHz luminosity of 100 /iJy kpc^). 

Figure [3] shows the fitted and measured positions for PSR J0108-1431. As shown in 
Table El the measured distance of 240^0^"^ pc is just consistent with the NE2001 distance, but 
inconsistent at 95% confidence with the earlier TC93 distance. It confirms that PSR J0108- 
1431 is more distant than originally suspected, and while its apparent luminosity is still low, 
it is no longer remarkably so. In fact, 5 other pulsars (PSR J0006+1834, PSR J1829+2456, 
PSR J1918+1541, PSR J2015+2524, and PSR J2307+2225) have lower apparent 400 MHz 
luminosities than the revised value of 510 yuJy kpc^ for PSR J0108-1431. 

These observations represent the first measurements of the proper motion of PSR J0108- 
1431, as well as the first distance measurement. Knowing the distance and transverse velocity 
of PSR J0108-1431 makes it possible to calculate the Shklovskii correction to its period 
derivative P, which is given by: 



Shk 

P 



where D is the pulsar distance, /i is proper motion and c is the speed of light. For PSR 
J0108-1431, the Shklovskii correction to the observed P is 1.37 x 10^^^, or 18% of the 



observed value of 7.704 x 10~^^ (iHobbs et al.ll2004l ). This correction to P allows a more 
accurate estimate of the pulsar's spin-down luminosity E = Att^IP/P^ and characteristic 
age Tc = P/ {{n — 1)P), where / is the neutron star moment of inertia (estimated to be lO'^^ 
g cm^) and n is the braking index, equal to 3 for magnetic dipole braking in a vacuum. For 
PSR J0108-1431, the value of E is revised from 5.8 x 10 
and the value of Tc is revised from 166 x 10^ years to 200 x 10^ years. 



30 gj,g g 1 to to 4.8 X 10^° erg s ^, 



PSR J0108-1431 has recently been detected in x-rays (jPavlov et al.l 120091) . prompting 
a re- examination of archival optical data and a subsequent optical detection (jMignani et al. 
20081 ) . The VLBI position and proper motion presented here confirms that the optical 
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detection of iMignani et al.l (120081 ) is coincident with the pulsar position at the time of the 
observation. The x-ray emission is consistent with a non-thermal spectrum of presumed 
magnetospheric origin, while the optical detection was consistent with thermal emission 
from the bulk of the neutron star surface. For both the x-ray and optical detections, the 
revised distance measure imphes a greater energy output in these wavebands, although the 
low precision of the parallax measurement means the luminosities remain relatively poorly 
constrained. 

For the x-ray waveband, the combination of revised distance measure and spin-down lu- 
minosity implies that PSR J0108-1431 actually convert s 1.5~^nj% of i t s spiii -down luminosity 



into x-rays, rather than the value of 0.4% reported in lPavlov et al.l (120091 ). This represents 
the highest known x-ray conversion efficiency of any pulsar - the distance and hence x- 
ray conversion efficiency of the previously most efficient pulsar (PSR J0630-2834) has been 
revised dramatically downwards by the measurements presented here. PSR J0630-2834 is 
discussed further in §4.31 

In the optical waveband, the upwardly revised luminosity requires either the addition 
of a large non-thermal component to the emission, or if the emission is solely thermal, 
a significant increase in the neutron star surface temperature. Any non-thermal optical 
emission would be generated at a much higher efficiency than other comparable old neutron 
stars, althoug h the ratio of o p tical efficiency to x-ray efficiency would be consistent with 
that found bv Izharikov et all J2004h . Unlike PSR B0950+08, the next oldest pulsar with 
multiwavelength data, the slope of a putati ve power-law compo nent could not be consistent 
with that inferred from x-ray observations (IMignani et al.ll2008l ). However, single power law 
fits are not ubiquitous amongst younger pulsars, and the contribution of non-thermal optical 
emission to the spectrum of PSR J0108-1431 cannot be discounted. 

If negligible magnetospheric optical emission is assumed, the implied surface tempera- 
ture for PSR J0108-1431 far exceeds that predicted for a neutron st ar of age 200 x 10^ yr 
(Shklovskii-corrected). Scaling the surface temperature obtained by IMignani et al.l (120081 ) 
by the increased optical luminosity at the best-fit distance of 240 pc yields a surface tem- 
perature of (2.9 ± 0.5) X 10^ K; taking the 95% confidenc e lower limit for distance (143 pc) 
and the lower temperature limit given in IMignani et al.l (|2008[ ) yields a minimum surface 
temperature of 8.5 x 10*^ K. According to standard neutron star cooling models, the maxi- 
mum temperature of an isolated neutron star of ag e 200 x 10^ yr should be 3 x 10'* K if no 
reheating occurs (see Figure 3 of lSchaab et al.lll999l ). To obtain a temperature above the the 
95% lower limit, the neutron star could be no older than 10 x 10^ yr - an unlikely prospect, 
as it would require the pulsar's b irth spin period to be at least 7 70 ms, much higher than 
typically assumed values (see e.g. iFaucher-Giguere &: Kaspill2006l ). Improved optical obser- 
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vations are needed to confirm or deny the presence of magnetosplieric optical emission, but 
tliese results are supportive of tlie presence of significant relieating for old, isolated neutron 
stars. 



4.2. PSR J0437-4715 



The results obtained for PSR J0437-4715 have already been discussed in iDeller et al. 



( 120081 ). and a brief summary is presented here. The fit to the observed pulsar positions is 
shown in Figure HI The pulsar distance of 156.3 ±1.3 pc which is derived from our VLBI 
astrometry is the most accurate pulsar distance (in both absolute and fractional distance) 
obtained for any pulsar to date. The comparison of this distance with the "kinematic dis- 



tance" derived by lVerbiest et al.l (120081 ) from the measurement of the rate of change of binary 



period yields a limit on the maximum appare nt acceleration of P SR J0437-4715 along the 



line of sight to the solar system. As shown in iDeller et al.l (120081 ). this can in turn be used 



to place limits on the maximum rate of change of Newton's gravitational constant G with 
time {G/G = (—5 ± 26) x 10^^'^ yr~^, 95% confidence), as well as the maximum mass of any 
undetected massive planets orbiting the solar system (no Jupiter-mass planets within 226 
AU of the solar system barycentre, 95% confidence over 50% of the sky), and the maximum 
amplitude of the characteristic strain spectrum of the stochastic gravitational wave back- 
ground (an amplitude of 1.1 x 10~^^ at a period of one year produces inconsistent distances 
in 95% of simulations). 

The covariance between proper motion and parallax for PSR J0437-4715 was the sec- 
ond highest of the pulsars in our sample after PSR J2048-1616. Since considerable time has 
passed since the end of our observing program, the addition of several new epochs could con- 
siderably improve the proper motion and parallax accuracy for PSR J0437-4715. Through 
simulations, we have estimated that three additional epochs over a year would reduce the 
parallax uncertainty for PSR J0437-4715 by a further factor of two. This extension to our 
observing program is the subject of a currently active LBA proposal. 



4.3. PSR J0630-2834 



PSR J0630-2834 is a mi ddle-aged pulsar which has bee n observed extensively as part of 



scintillation studies (see e.g. ICordeslll986l : iBhat et al.l 119991) . and wou l d be undistinguished 



if not for its remarkable apparent luminosity in x-rays. iBecker et al.l (120051 ) observed PSR 



J0630-2834 with XMM-Newton and found an x-ray luminosity of 8.4 x 10 erg s , based 
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on the NE2001 distance of 1.45 kpc. This imphed that the pulsar was converting ~ 16% 
of its spin-down luminosity into x-rays - an order of magnitude more than any other old 
or middle-aged pulsar. They suggest that the most likely explanation is that the distance 
to the pulsar is over-estimated, a theory which has been proven correct by these VLBI 
observations. 

The measured and fitted positions of PSR J0630-2834 are shown in Figure O The 
measured distance of 332I4Q pc, as shown in Table [21 means that the actual x-ray conversion 
efficiency is a much less surprising 

O.Stoi %■ It is also within the la error bars of the efficiency 
derived for PSR J0108-1431 in §4.11 above. The best-fit distance of 332 pc is less than a 
quarter of the DM-derived distance in the NE2001 model, and a factor of seven smaller than 
the distance predicted by the TC93 model. This large discrepancy is discussed further in 

The proper motion of PSR J0630-2834 has previously been measured using the VLA 



by iBrisken et al.l (j2003af ). and the measured values from this work (/Zq, = —46.3 ± 1.0 mas 
yr~^. Us = 21.3 ± 0.5 mas yr~^) are consistent with, but marginally more precise, than the 
VLA results (/i^ = — 44.6±0.9mas yr~^, fis = 19.5±2.2mas yr^^). The resultant Shklovskii 
correction is less than 0.04% of the observed period derivative (7.12 x 10~^^). 



4.4. PSR J0737-3039 



The results obtained for PSR J0737-3039 have already been discussed in iDeller et al. 



(j2009al ). and a brief summary is presented here. The fit to the observed pulsar positions 



is shown in Figure O The measured distance of 1150^*^^20 pc inconsistent with previous 
DM-based estimates of 570 pc (TC93) and 480 pc (NE20 01), and is much mo re significant 
than the marginal timing parallax detection of 333^*^^33 pc ( Kramer et al.ll2006l ). 



The kinematic information obtained from these measurements can be used to accu- 
rately calculate the Shklovskii term and other kinematic contributions to the binary period 
derivative of PSR J0737-3039A/B for the first time, allowing an estimation of the precision 
to which GR can be tested. The measurements reported here reduce the uncertainty on 
kinematic terms in the observed binary period derivative to 0.01% of the GR contribution, 
implying that with increasing accuracy of timing observations, tests of GR can be made 
in the strong-field regime to the 0.01% level. The kinematic information also implies that 
the transverse velocity of PSR J0737-3039A/B in the local standard of rest (LSR) is only 
913 km s~^. Given that the system must have survived two supernovae undisrupted, this low 
velocity is surprising, and lends support to models of formation of the system in which the 
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neutron stars receive low or no kicks at birth (e.g. iPiran fc Shavivll2005l ) 



Finally, the revised distance measure helps to clarify the origin of the bulk of the x-ray 
emission in the PSR J0737-3039A/B system. A magnetospheric origin for the x-radiation 
(from pulsar A) had been considered a strong possibility, but the hydrogen column density 
calculated from fits to the x-ray data had appeared to be unreasonably h igh, given the 



presu med location of PSR J0737-3039A/B in the Gum nebula HII region (jPossenti et al. 



20081 ). The revised distance measure places the pulsar beyond the Gum nebula and makes 



the calculated hydrogen column density consistent with expectations. 



4.5. PSR J1559-4438 



The results obtained for PSR J1559-4438 have already been discussed in iDeller et al. 



( l2009bl ). and a brief summary is presented here. The fit to the observed pulsar positions 
is shown in Figure El At a distance of 26OOI45Q P*^' PSR J1559-4438 is one of the most 
distant pulsars with a parallax measur ed using VLBI. The rn easured distance is consistent 
with the NE2001 prediction (2350 pc; ICordes fc Lazid l2002l ). which differed considerably 
from the earlier TC93 distance estimate of 1580 pc. It is also c onsistent with th e lowe r 
distance estimate of 2.0 ±0.5 kpc made using HI line absorptions by iKoribalski et al.l (119951 ). 
The measured values of proper motion (/i^ = 1.52 ± 0.14 mas yr~^, /^^ = 13.15 ± 0.05 mas 
yr~^) are consistent with the VLA observations of iFomalont et al.l (119971 ). who measured 
/Iq = 1 ± 6 mas yr~^, = 14 ± 11 mas yr~^. 

With an accurate proper motion now calculated, the position angle of the proper motion 
for PSR J1559-4438 can be compared to the position angle of the emissi on polarisation. 



which tests the alignment of the rotation and velocity vector, as described by I Johnston et al. 



( 120071 ). In the case of PSR J1559-4438, the position angle of the proper motion (6.6 ± 0.6 °) 
is not significantly aligned para llel with or perpend icular to the polarisation position angle 
(71 ± 3 °). However, as noted i nlDeller et al.l (2009b ). the polarisation profile of PSR J1559- 



4438 (as shown in Figure 5 of I Johnston et al.ll2007l ) is complicated, and it is possible that 
the original determination of the magnetic field orientation from the polarisation position 
angle was incorrect. 



4.6. PSR J2048-1616 



PSR J2048-1616 was the second "technique check" source included in the observing 
program, after PSR J1559-4438. It is less bright than PSR J1559-4438, but was predicted 



- 11 - 



to be somewhat closer (640 pc in the TC93 model, 560 pc in the NE2001 model). As 
shown in Tabled whilst a parallax was measured for PSR J2048-1616, it was not significant 
(1.9cr). Nevertheless, the best-fit distance of 580 pc is consistent with the DM-based distance 
estimates, and an accurate measurement of proper motion was made. Figure [H] shows the 
fitted and measured positions of PSR J2048-1616. 

It is apparent from Figure [8] that the large position error of the third epoch (MJD 
54182) is the primary reason for the poorly constrained fit. Equipment failure at the Mopra 
telescope and limited time on the Tidbinbilla telescope during this epoch reduced the number 
of telescopes on-source during observations of PSR J2048-1616 to three, and thus the large 
errors are unsurprising. The covariance between proper motion and parallax was the largest 
for PSR J2048-1616, and hence one or more additional position measurements, preferably 
close to the appropriate parallax extrema, could have significantly improved the parallax 
determination. 

A considerably more accurate parallax for PSR J2048-1616 has recently been reported 
by lChatterjee et al.l (120091 ) . who observed PSR J2048- 1616 along with 13 other pulsars using 



the VLBA. The parallax of 1.05+°;°^ mas obtained by IChatterjee et al.l (120091 ) is consistent 



with our results (1.71±0.91 mas), but the proper m otion values differ sign ificantly (/Xq, 



113.16 ± 0.02 mas ji'^, fig = -4.60+0:1 mas yr'^ for IChatteriee et al.ll2009l : /i^ = 114.24 ± 



0.52 mas yr~^, fis = —4.03 ± 0.24 mas yr~^ from our observations). Combining the epoch 
position fits from both observing campaigns and allowing a constant offset between the LBA 
and VLBA epochs yields a fit which is consistent with the VLBA-only values, with similar 
errors. This is unsurprising since the VLBA epochs, with much lower errors, dominate the fit. 
Refitting the LBA-only data after fixing the parallax to the VLBA-derived value of 1.05+o o2 
mas removes the impact of the covariance with parallax and yields fia = 113.94 ± 0.32 mas 
yr^^ and = —4.30 ± 0.10 mas yr~^; closer to the VLBA-only values, but still differing at 
the 2a level in right ascension. 

Whilst both observing campaigns utilized t he same reference sour ce (J2047-1639), they 
sampled different time periods (2002 - 2005 for IChatteriee et al.l[2009l . 2006 - 2008 for these 
observations) and the VLBA observations utilized an in-beam calibrator. The primary 
calibrator is slightly resolved, as shown in Figure [H and thus imperfect modeling of its 
(possibly time varying) structure is a plausible explanation for the proper motion discrepancy. 
The problem of modeling the structure of resolved calibrator sources is particularly acute 
for the LBA, due to the limited uv coverage. 
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4.7. PSR J2144-3933 



^PSR J2144-3933 was discovered in the Parkes Southern Pulsar S urvey (IManchester et al 

1996 ). and was initially misidentified with a period of 2.83 seconds. Young et al. ( 19991 ) re- 
ported that the true period is in fact 8.5 seconds, making PSR J2144-3933 the longest 
period radio pulsar k nown. PSR J2144-3933 hes below the traditionally assumed pulsar 
"death line" (see e.g. IChen fc RudermanI 119931 ). and hence its true luminosity is extremely 
important for models of pulsar emission and evolution. The discovery of PSR J2144-3933 
prompted alternative models of pulsar emission in w hich long period p ulsars such as PSR 
J2144-3933 could remain luminous in the radio (e.g. Izhang et al.lboooh . PSR J2144-3933 
possesses a steeper spectral index than average (—2.4), and while several other pulsars are 
less luminous at 400 MHz, PSR J2144-3933 is the least luminous pulsar known at 1400 MHz 
(24 /iJy kpc^ at assumed distance of 180 pc; average 1400 MHz flux density of 0.75 mJy 
calculated from archival Parkes observations). 

The fitted and measured positions of PSR J2144-3933 are shown in Figure O As 
shown in Table [3], a highly significant parallax was detected, corresponding to a distance of 
165l}4pc. Given the generally assumed errors on DM distances, this is consistent with the 
TC93 value (180 pc), but not the NE2001 value (264 pc). This confirms that the appar- 
ent radio luminosity of PSR J2144-3933 is extremely low - 15% lower than the previously 
assumed value of 24 fiJj kpc^. 

The proper motion measurement allows for the first time a calculation of the Shklovskii 
correction for the period derivative of PSR J2144-3933 - the Shklovskii contribution to P is 
9.4 X 10^^'', or approximately 19% of the observed P of 4.96 x 10^^^. The true spin-down 
luminosity of PSR J2144-3933 is thus further reduced from the assumed value of 3.2 x 10^^ erg 
s~^ to 2.6 X 10^^ erg s~^, the smallest known spin-down lur ninosity of any pulsa r. PSR J0343- 
3000, discovered in the Parkes High-Latitude survey by iBurgay et al.l (120061 ) . has the next 
lowest spin-down luminosity, which is a factor of 5 higher than the revised value for PSR 
J2144-3933. In addition, this revision to the P value for PSR J2144-3933 places it even 
further past the assumed pulsar "death line" . 



As noted by lYoung et al 



(119991 ). the low luminosity and narrow beaming fraction 
(~ 0.01) of PSR J2144-3933 imply that many such objects exist in the Galaxy. These 
measurements confirm that the true distributions of pulsar spin-down luminosities and ra- 
dio luminosities do extend down to very low levels, which is important for the generation 
of synthetic pulsar catalogues and predictions of discovery rates for instruments such as the 
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Square Kilometre Array (SKA@) and LOFAR0. Assuming the steep spectral index of PSR 
J2144-3933 is typical for similar long-period, low-luminosity objects, and that the Galactic 
population of such objects can be estima ted from PSR J214 4-3933 alone, more than 100,000 
similar pulsars may inhabit the Galaxy (lYoung et al.lll999l ). 



4.8. PSR J2145-0750 

PSR J2145-0750 is a binary pulsar with s pin p eriod 16.05 ms and orbital p eriod 6.84 



days, which was discovered by iBailes et al.l (jl994l ). Optical observations by iBell et al. 



( 119951 ) tentatively identified the expected white dwarf companion to PSR J2145-0750, which 
was confirmed bv iLundgren et al.l (119961 ). Timing astrometry using a 10 year dataset by 
Lohmer et al.l (120041 ) measured a parallax of 2.0 ± 0.6 mas, and a proper moti on of 14.1 ± 
0.4 ma s yr~^, corresponding to a transverse velocity of 33 ± 9 km s~^. However, iHotan et al. 



(120061) reported no detection of parallax and an upper limit of 0.9 mas, at 95% confidence, 
albeit with a shorter (2.5 year) dataset, in addition to a proper motion measurement of 
13.5 ± 6.0 mas yr ~^. A previous atte mpt has been made to measure the VLBI parallax 
of this system by iBrisken et al.l (|2002| ) using the VLBA, but the pulsar was not detected 
and dropped from the observing program. Accordingly, the distance to this system is still 
controversial and confirmation of a previous result is keenly sought. 

From the four observations made of PSR J2145-0750 during this observing program, 
significant detections were mad e on only two occasions. As was the case for the VLBA 
program of iBrisken et al.l (120021 ). it is believed that refractive scintillation is the major cause 
of the non-detections, as the pulsar's flux varies significantly over long timescales. Although 
measurement of the parallax of PSR J2145-0750 was thus not possible, by holding parallax 
fixed at zero mas it was possible to measure a proper motion of Ha = —15.4 ± 2.1 mas yr~^. 
Us = —7.7 ± 0.80 mas yr~^. Varying the parallax between zero and two mas results in a 
change in proper motion of less than 0.1 mas yr~^, much smaller than the formal errors. As 
the proper motion and position have been derived from only two measurements, there are no 
degrees of freedom to the fit and no estimate of systematic errors, and so the result should 
be treated with some caution. The fitted and measured positions of PSR J2145-0750 are 
shown in Figure [TOl 

The observed total proper motion of 17.2 ± 2.2 mas yr~^ is consistent with the value of 



^http: / /www. skatelescope.org 
^http: / /www. lofar.org/ 
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13.5 ±6.0 mas yr ^ obtained by lHotan et al.l ( 20061). but only cons istent at the 2o level with 
the value of 14.1 ± 0.4 mas yr~^ obtained by Lohmer et al. (2004]§. As these VLBI results 



are derived from only two measurements, and there is no way of quantifying the systematic 
error, it is likely that the lower value of proper motion derived from timing (~14 mas yr~^) 
is more accurate for this pulsar. 



4.9. Comparison to DM distance models 

Galactic electron density distribution models are extremely important for all fields of 
pulsar research, since their prediction of pulsar distances based on DM has the potential 
to bias estimates of luminosity, velocity, and various timing terms for the vast majority of 
pulsars without an independent distance constraint. Hence, continual improvement of these 
models, and characterisation of errors in the existing models (and the impact of these errors) 
is a crucial task. 

The seven pulsar parallaxes presented here make a substantial addition to the 44 pulsars 
with published parallaxes, and hence a review of the accuracies of the TC93 and NE2001 
models is timely. For the comparison presented below, only pulsars with a parallax more 
significant than 2cr are considered, and where multiple parallax measurements exist, only the 
most accurate is used. This leads to a sample size of 41 (33 from VLBI and 8 from timing). It 
is appropriate to note that large distance errors are likely to be over-represented in existing 
astrometric results. One reason is the selection effect of anomalous pulsars being chosen 
for study, which was certainly the case for the PSR J0630-2834 system. Another potential 
factor is that most astrometric distance determinations to date have been for relatively 
nearby pulsars, where there is less opportunity for underdensities and overdensities in the 
ISM along the line of sight to cancel, and hence it is plausible that distance models are on 
average more reliable for more distant pulsars. 

Figure [TT] shows a histogram of the observed distance errors for the TC93 model (where 
the best-fit parallax distances have been assumed to be the correct value), with the error ex- 
pressed in dB, and binned in 1 dB increments. The error in dB is defined as 10 logm '""ff^ distance 

^ ' ^J^u parallax distance ' 

and hence positive values represent overestimates, and negative values represent underesti- 
mates. Figure [12] repeats this plot for the NE2001 model. Figure [TT] shows a clear systematic 
offset in the median error of the TC93 model towards underestimated distances, but the 



^Since PSR J2145-0750 lies close to the ecliptic, timing observations obtain large covariances between 
proper motion in right ascension and declination, and so the comparison of total proper motion is more 
useful. 
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largest errors are seen when distances are overestimated. Figure [12] shows that systematic 
bias has been largely removed by the NE2001 model, but the distribution of errors still can- 
not be well approximated by a single Gaussian model, due to the long wings towards large 
errors. For both the TC93 and NE2001 models, the largest errors are seen when the distance 
is overestimated. The effect of these high-error "outliers" on population studies of neutron 
stars is considered further below. 

Estimating the true form of the error distribution shown in Figure [12] is difficult due to 
the small number of samples. However, the single Gaussian model with standard deviation 
0.8 dB plotted as a solid line in Figure [12] shows that the common assumption of 20% errors 
for DM-based distance estimates is not realistic. In the discussion below, this model (in 
which 67% of errors are less than a factor of 1.2) is referred to as the traditional error model. 

Rather than attempting to discern the true error distribution, we have modeled the 
distribution by simply summing the probability distribution function of the distance error for 
each pulsar, which accounts for the possible variation of the true distance from the parallax 
value. The total probability distribution function was binned in 0.001 dB increments and 
the result smoothed with a 1 dB Banning window. This is plotted as a dashed line in 
Figure [12], and is hereafter referred to as the binned error model. The inclusion of parallax 
uncertainties acts to broaden the distribution slightly, and somewhat preferentially towards 
more negative errors. This is caused by a uniform parallax error encompassing a non-uniform 
distance space, with a greater range of distances if the parallax was overestimated t han 



underestimated (a manifestation of the so-called Lutz-Kelker bias; iLutz fc Kelkerlll973l ) 



Three pulsars with distance errors >5.5 dB dominate the tail of the NE2001 error 
distribution with overestimate d distances. They are PSR J0613-0200 (NE2001 distance 



1700 pc, parallax 2.1 ±0.6 mas. lHotan et al.ll2006f ). PSR J0630-2834 (NE2001 distance 1450 



pc, parallax 3.0 ± 0.4 m as, this work) and PS R B1541+09 (NE2001 distance 35000 pc 



parallax 0.13 ± 0.02 mas, IChatteriee et al.ll2009l ). The NE2001 distance of 35 kpc for PSR 



B1541+09 is somewhat arbitrary, as the NE2001 model does not possess sufficient electrons 
along the line of sight to PSR B1541+09 to account for its DM, and hence it is placed beyond 
the edge of the model. This type of shortcoming is rare and can be easily identified, and 
so we have excluded PSR B1541+09 from the analysis that follows. Figure [13] shows the 
revised error histogram and binned model. Despite the exclusion of PSR B1541+09, PSR 
J0630-2834 is still an example of the previously mentioned selection effect bias, and so (as 
already noted) large errors are likely to be somewhat over-represented in the subsequent 
analysis. Thus, this model can safely be considered to be a "worst-case" representation of 
DM distance reliability. 
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4.10. The pulsar velocity distribution 

Although is has been been shown that there is no evidence for systematic bias in the 
NE2001 distance model, errors exceeding 6 dB still exist for individual pulsars. Such errors 
can impact upon population analysis of pulsars, artificially creating tails of extremely high or 
low values of distance-dependent parameters such as luminosity and velocity. This can affect 
estimates of the mean values of these quantities, as well as confusing attempts to explain 
the underlying physics by generating false outliers. As an example, the case of neutron star 
transverse velocities is considered below. 



Hobbs et al.l (120051 ) examine a large sample of pulsar proper motions and conclude that 
the distribution of 3D space velocities of young pulsars (characteristic age < 3 Myr) is well fit 
by a Maxwellian distribution with mean 431 km s~^ and one-dimensional rms 265 km s^^. 
This velocity distribution is used as the starting point for the simulations that follow. A 
Monte Carlo simulation was performed, creating a synthetic pulsar catalogue of ten million 
pulsars, where each pulsar's actual velocity was drawn from the distribution described above. 
These are referred to as the unperturbed velocities. No proper motion measurement error 
was assigned to the unperturbed velocities; more distant pulsars (which possess on average a 
lower angular proper motion) would be more affected by the presence of a roughly constant 
proper motion measurement error, which is another small but unmodeled bias which would 
act to further broaden the observed velocity distribution. The velocities of each pulsar were 
perturbed according to the binned and the traditional distance error functions described 
above, and the observed 2D velocity was recorded for each case. The observed 2D velocity 
distribution for the binned error function and the traditional error function, along with the 
unperturbed 2D velocity distribution, is plotted in Figure HH 

It is immediately apparent that the broadening of the observed velocity distribution is 
considerably greater for the binned error model, compared to the traditional error model, 
and is particularly pronounced at the high end of the velocity distribution. The mean 
unperturbed 2D velocity is 332 km s~^, which is increased by 2% to 338 km s~^ in the 
traditional error model, and by 5% to 349 km s~^ in the binned error model. 

More important than the mean effect on pulsar veloci t ies, h owever, is the effect on the 



high end of the pulsar velocity distribution. iHobbs et al.l (120051 ) note that distance model 



inaccuracies could be responsible for the sparse tail of very high velocity pulsars, but conclude 
that with the small sample size available, the observed high-velocity pulsars are consistent 
with the tail of a continuous velocity distribution. In considering the impact of distance 
model errors, it is useful to examine the frequency of occurrence of pulsars with observed 
transverse speeds greater than 1000 km s~^ in the synthetic catalogue. These pulsars will 
be described as very high velocity (VHV) pulsars. In the synthetic catalogue, only 0.1% of 
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pulsars fall into this category in the unperturbed velocity distribution, but the occurrence 
rises to 0.4% when the traditional error model is applied, and 3.9% for the binned error 
model. By way of comparison, the ATNF pulsar catalogue contains 216 pulsars with listed 
transverse velocities - of these, seven exceed 1000 km s~^, approximately 3% of the sample. 

Thus, distance model errors can add to the high-velocity tail observed in the pulsar 
velocity distribution, an effect which is significantly enhanced when the distance errors "out- 
liers" are accounted for. These results suggest that the majority of VHV pulsars may in 
fact simply be pulsars with incorrect distance estimates. As already stated, the binned error 
model used here is likely to overestimate the true frequency of large errors, and hence the true 
impact of DM distance errors is likely to lie somewhere in between the binned error model 
and the traditional error model. Additio nally, there is clear evid ence that at least some VHV 
pulsars are present in the Galaxy, with IChatterjee et al.l (120051 ) using VLBA astrometry to 
show that the transverse velocity of PSR B1508-I-55 exceeds 1000 km s~^, and observations 
of PSR B2224+65 showing scintillat ion, proper r notion and bow shock results consisten t 
with a velocity close to 1000 km s~^ (ICordeslll986l : [Harrison et al.lll993l : ICordes et al.lll993l ). 
Nonetheless, the binned error model is useful in illustrating that the largest impact of occa- 
sional, rare distance errors is at the high end of the pulsar velocity distribution. A larger, 
unbiased sample of pulsar distance measurements would enable a more accurate estimation 
of the true distance error function for NE2001. This would in turn enable the effects of dis- 
tance errors to be "deconvolved" from the measured pulsar velocity distribution with some 
confidence, and would lay the framework for a further improved Galactic electron density 
distribution model to supersede NE2001. 



5. Conclusion 



Using the Long Baseline Array, we have measured parallaxes for seven pulsars in the 
sparsely sampled southern sky, a four-fold increase on the previous number of southern 
VLBI parallaxes. Based on the results of this survey, we conclude that parallaxes could 
be measured for at least 30 additional southern pulsars using the LBA in its present form, 
and future enhancement to LBA sensitivity can be expected to allow parallax measurements 
for an increasing sample of fainter pulsars. The ensemble of independent distance measures 
obtained from this program will allow significant improvements to Galactic electron models at 
these latitudes. Two of the measured distances (PSR J0630-2834 and PSR J0737-3039A/B) 
differed by more than a factor of two from their DM estimates, reinforcing that DM-based 
distance estimates for individual pulsars should be treated cautiously in the absence of any 
other distance indicators. Finally, it was shown that DM distance errors could be inflating 
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estimates of both the mean pulsar velocity and the number of very high velocity pulsars, 
although the current small and biased sample of independent distance measures makes this 
effect difficult to quantify. 

The authors would like to thank S. Chatterjee for providing the fitted positions of PSR 
J2048-1616 from VLBA observations. The Long Baseline Array is part of the Australia 
Telescope which is funded by the Commonwealth of Australia for operation as a National 
Facility managed by CSIRO. This research made use of the ATNF Pulsar Catalogue and 
the NASA/IPAC Extragalactic Database (NED). NED is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 
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Table 1. Target pulsars 



Pulsar 


TC93 distance 


NE2001 distanee 


1600 MHz 


Pulsar 


Equivalent gated 


R.efercncc 


Calibrator / target 


name 


(pc) 


(pc) 


flux (mjy) 


gating gain 


1600 MHz flux (mJy) 


source 


separation (dcg) 


PSR J0108-1431 


130 


180 


0.8" 


5.1 


4 


JOlll-1317 


1.5 


PSR J0437-4715 


140 


190 


140 


6.25 


875 


J0439-4522 


1.9 


PSR J0630-2834 


2150 


1450 


23 


3.5 


81 


J0628-2805 


0.7 


PSR J0737-3039A/B 


570 


570 


1.6 


2.5 


4 


J0738-3025 


0.4 


PSR J1559-4438 


1600 


2350 


40 


3.6 


144 


J1604-4441 


0.9 


PSR J2048-1616 


640 


560 


13 


3.6 


47 


J2047-1639 


0.5 


PSR ,72144-3933 


180 


264 


0.8" 


10 


8 


J2141-3729 


2.1 


PSR ,72145-0750 


500 


570 


8" 


4.3 


34 


,72142-0437 


3.3 



^Pulsar suffers heavily from long timescale scintillation, so individual epochs vary considerably from the average value shown. 



Table 2. Reference source positions 



Source Right ascension (J2000) Declination (J2000) Reference 



JOlll- 


-1317 


Ol'"!!" 


'56.857962' 


-13 


°17'01. 19702" 


VCS5: Kovalev et al. f2007) 


J0439- 


-4522 


04''39" 


'00.854714= 


-49 


°22'22. 56260" 


ICRF: Ma et al. ri9981 


J0628- 


-2805 


06''28" 


'43.279032= 


-28 


°05'19. 38364" 


VCS5; Kovalev et al. nOOT) 


B0736 


-303 


07'>38" 


'19.78785= 


-30 


°25'04.8360" 


None: see Deller et al. f2009al 


J1604- 


-4441 


16'' 04" 


'31.020692= 


-44 


°41'31. 97421" 


ICRF; Ma et al. ri998) 


J2047- 


-1639 


20''47" 


'19.667018= 


-16 


°39'05. 84249" 


VCSl: Beaslev et al. f20021 


J2141- 


-3729 


2l''41" 


'52.448964= 


-37 


°29'12. 99259" 


VCS3; Potrov et al. f2005) 


J2142- 


-0437 


2l'>42" 


'36.901688= 


-04 


°37'43. 51275" 


VCSl; Beaslev et al. (20021 



Table 3. Astrometric fits for all target pulsars 



PSR ,70108-1431 



PSR J0437-4715 



PSR J0630-2834 



PSR .10737-3039 



Right ascension (J2000) 




01:08:08.347016 ± 0.000088 


04:37:15.883250 ± 0.000003 


06:30:49.404393 ± 0.000043 


07:37:51.248419 ± 


0. 


000026 


Declination (J2()()()) 




-14:31:50.187139 ± 0.001069 


-47:15:09.031863 ± 0.000037 


-28:34:42.778813 ± 0.000372 


-30:39:40.714310 ± 


0. 


000099 


Proper motion in right ascension 


(fj^a ', mas yr "'" ) 


75.05 ± 2.26 


121.679 ± 0.05 


-46.30 ± 0.99 


-3.82 ± 


0. 


62 


Proper motion in declination (/Li^ 


; mas yr""*^) 


-152.54 ± 1.65 


-71.820 ± 0.09 


21.26 ± 0.52 


2.13 ± 


0. 


23 


Parallax (tt; mas) 




4.170 ± 1.421 


6.396 ± 0.054 


3.009 ± 0.409 


0.872 ± 


0. 


143 


Distance (pc) 






156.3+}:'j 


332j:^^ 








Transverse velocity {vf, km s""'") 




194±lf 


104.7;;;;^ 


80115 


24t^ 







Reference epoeh (MJD) 
Visibility weigliting 
Image weighting 

Average epocii mean fit error (mas) 
Average intra— epoclr systematic error (mas) 
Average inter— epoch systematic error (mas) 
Average single— epoch SNR 
Covariance of tt and fj-^t 
Covariance of tt and 



54100.0 
Sensitivity 
Natural 
1.232 
2.477 
4.310 
8 

0.27 
0.00 



54100.0 
Equal 
Uniform 
0.059 
0.068 
0.103 
21 
-0.39 
0.60 



,54100.0 
Equal 
Uniform 
0.765 
0.839 
1.205 
15 
-0.01 
-0.10 



54100.0 
Sensitivity 
Natural 
0.747 
0.939 
0.0 
17 
-0.01 
-0.06 



Parameter 



PSR, .11559-4438 



PSR ,12048-1616 



PSR ,12144-3933 



PSR .7214,5-0750* 



to 

CO 



Right ascension (72000) 
Declination (J2000) 

Proper motion in right ascension {^ai mas yr~^) 
Proper motion in declination (fJ-g; mas yr~^) 
Parallax (tt; mas) 
Distance (pc) 

Transverse velocity (vt', km s""*^) 
Reference epoch (MJD) 
Visibility weighting 
Image weighting 

Average epoch mean fit error (mas) 
Average intra— epoch systematic error (mas) 
Average inter— epoch systematic error (mas) 
Average single-epoch SNR 
Covariance of 7r and ^ct 
Covariance of tt and ^§ 



n 4-690 

-450 
0-F44 



15:59:41.526126 ± 0.000008 
-44:38:46.901778 ± 0.000035 
1.52 ± 0.14 
13.15 ± 0.05 
0.384 ± 0.081 
2600l 
1631 
54100.0 
Equal 
Uniform 
0.242 
0.259 
0.055 
50 
-0.43 
0.18 



20:48:35.640637 ± 0.000040 
-16:16:44.553501 ± 0.000147 
114.24 ± 0.52 
-4.03 ± 0.24 
1.712 ± 0.909 
_„„+660 
„, -4-362 

■31'_111 

54100.0 
Equal 
Uniform 
0.617 
1.282 
0.106 
23 
0.78 
0.64 



12:44:12.060404 ± 0.000045 
-39:33:66.885041 ± 0.000316 

-57.89 ± 0.88 
-155.90 ± 0.54 
6.051 ± 0.660 



5+17 



165-14 

i3oi;-, 

54100.0 
Equal 
Natural 
1.026 
1.460 
0.876 
10 
0.23 
0.16 



21:45:50.461901 ± 0.000098 
-07:50:18.462388 ± 0.000668 
-16.43 ± 2.07 
-7.67 ± 0.81 



54100.0 
Sensitivity 
Natural 
2.136 
0.0 
0.0 



^ Based on two detections, with parallax fixed at 0. Variation of the parallax value between and 2 mas results in less than 100 juas yr ^ difference in derived proper motion. 
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Right Ascension (mas) Right Ascension (mos) 

Fig. 1. — Images of the phase reference sources used. In aU instances, the image contours 
increase by factors of two from a minimum value expressed as a percentage of peak image 
flux. Runnning left to right, top to bottom, the calibrator name, peak image flux (mJy) 
and minimum contour are: JOlll-1431, 188, 0.5%; J0439-4522, 440, 0.1%; J0628-2805, 117, 
0.5%; B0736-303, 116, 1%. The image of B0736-303 has been shifted to place the image 
centre between the two components for clarity. 
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Fig. 2. — Images of the phase reference sources used. In all instances, the image contours 
increase by factors of two from a minimum value expressed as a percentage of peak image 
flux. Runnning left to right, top to bottom, the calibrator name, peak image flux (mJy) and 
minimum contour arc: J1604-4441, 868, 0.25%; J2047-1639, 336, 0.5%; J2141-3729, 320, 
0.25%; J2142-0437, 280, 1%. 
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Fig. 3. — Motion of PSR J0108-1431, with measured positions overlaid on the best fit. 
From top to bottom; motion in dechnation vs right ascension, motion in right ascension vs 
time with proper motion subtracted, and motion in dechnation vs time with proper motion 
subtracted. 
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Fig. 4. — Motion of PSR J0437-4715, with measured positions overlaid on the best fit. 
From top to bottom; motion in dechnation vs right ascension, motion in right ascension vs 
time with proper motion subtracted, and motion in dechnation vs time with proper motion 
subtracted. 
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Fig. 5. — Motion of PSR J0630-2834, with measured positions overlaid on the best fit. 
From top to bottom; motion in dechnation vs right ascension, motion in right ascension vs 
time with proper motion subtracted, and motion in dechnation vs time with proper motion 
subtracted. 
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Fig. 6. — Motion of PSR J0737-3039A/B, with measured positions overlaid on the best fit. 
From top to bottom; motion in dechnation vs right ascension, motion in right ascension vs 
time with proper motion subtracted, and motion in dechnation vs time with proper motion 
subtracted. 
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Fig. 7. — Motion of PSR J1559-4438, with measured positions overlaid on the best fit. 
From top to bottom; motion in dechnation vs right ascension, motion in right ascension vs 
time with proper motion subtracted, and motion in dechnation vs time with proper motion 
subtracted. 
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Fig. 8. — Motion of PSR J2048-1616, with measured positions overlaid on the best fit. 
From top to bottom; motion in dechnation vs right ascension, motion in right ascension vs 
time with proper motion subtracted, and motion in dechnation vs time with proper motion 
subtracted. 
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Fig. 9. — Motion of PSR J2144-3933, with measured positions overlaid on the best fit. 
From top to bottom; motion in dechnation vs right ascension, motion in right ascension vs 
time with proper motion subtracted, and motion in dechnation vs time with proper motion 
subtracted. 
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Fig. 10. — Motion of PSR J2145-0750 in right ascension plotted against motion in declina- 
tion, with measured positions overlaid on the best fit. Since only two significant measure- 
ments were made, the parallax of the pulsar was held fixed at zero, and the fit to position and 
proper motion has no free parameters. As the parallax was fixed at zero, the plots of motion 
in right ascension vs time and declination vs time which were shown for other pulsars are 
not shown here. Variation of the parallax between zero and two mas made an insignificant 
(< 100 //as yr~^) difference to fitted proper motion. 
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Fig. 11. — Histogram of TC93 errors for pulsars with measured parallaxes, with errors binned 
in 1 dB increments. Underestimates of distance arc more common in this model, but the 
largest errors are made when the distance is overestimated. The standard deviation of the 
errors is 2.8 dB. 
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Fig. 12. — Histogram of NE2001 errors for pulsars with measured parallaxes, with errors 
binned in 1 dB increments. Unlike the TC93 model, there is no discernable systematic 
errors, but the distribution is still clearly non-Gaussian, with a long tail of errors which can 
exceed 6 dB. As with TC93, the largest errors are seen when the distance is overestimated. 
The standard deviation of the errors is 2.3 dB, a small improvement on the 2.8 dB seen for 
TC93. The "binned" error model described in the text is shown as a dashed line, while the 
single Gaussian "traditional" model which approximates the 20% errors commonly assumed 
for DM distances is shown as a solid line. 
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Fig. 13. — Histogram of NE2001 errors for pulsars with measured parallaxes, with errors 
binned in 1 dB increments, after the exclusion of PSR B1541+09. The "binned" error model 
is shown as a dashed line, and the "traditional" error model is shown as a sohd line. 
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Fig. 14. — 2D velocities for a synthetic pulsar population, as observed with no distance error 
(solid line), the traditional distance error model (dashed) and the binned distance error 
model described in the text (dash-dot). The effect of the binned distance error model is to 
broaden the observed velocity distribution, particularly at large velocities. 



